Cation site preference in zeolite clinoptilolite: A density functional study

The structure of zeolite clinoptilolite and the coordination of extraframework cations (Ca, Ba, Na and K) by framework oxygen atoms are examined by density functional theory (DFT) using two methods: the ONIOM two-layer model (DFT/MM) and periodic DFT. Both ONIOM and periodic model calculations predict correctly the clinoptilolite structure and find configurations with Al ? Si substitution at the T2 site which interconnects the xz layers via its apical oxygen atom as the most stable ones, in agreement with experimental data. The configurations with Al? Si substitution at the T1 sites are favorable for migration of cations into the large channel A from the eight-member rings which are side rings of channel A and opened in channel C (c-rings); the process is energetically most favorable for the Na cations. K and Ba cations occupy sites in the vicinity of the c-rings: M1, M1a and M3, among which K prefers the M3 site inside the c-rings. The molecular electrostatic potential (MEP) maps reveal areas of increased nucleophilic properties inside the large channel A, to which cations can migrate upon heating, or incoming cations can be retained. The minimal interatomic distances between two equivalent cations, residing in channel A are 4.850 Å (Ca–Ca); 5.025 Å (Na–Na) and 4.305 Å (K–K). The ONIOM method should be preferred over periodic models for describing cations in the c-rings: the smaller c-parameter of the unit cell replicates the cations in adjacent rings along the [001] direction. 2013 Elsevier Inc. All rights reserved.